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Abstract 

The many body contributions to the nuclear anapole moment of 133 Cs, 205 T1, 
207 ' 209 PB, and 209 Bi from the core polarization are calculated in the random-phase 
approximation with the effective residual interaction. Strong reduction of a valence 
nucleon contribution was found provided by the core polarization effects. The con- 
tribution of the core particles to the anapole moment compensates this reduction to 
large extent keeping the magnitude of nuclear anapole moment close to its initial single 
particle value. 



1 Introduction 



The atomic parity non- conservation (PNC) effects dependent on nuclear spin are expected 
to be dominated by contact electromagnetic interaction of electrons with nuclear anapole 
moment (AM) |2[. The anapole is a new electro- magnetic moment arising in a system 
without center of inversion ||. It exists even in such a common object as a chiral molecule in 
a state with nonvanishing angular momentum |3J. The nuclear anapole moment is induced 
by PNC nuclear forces. 

In all calculations of the anapole moment M |5], |7|, an independent particle model has 
been used for the nucleus. In this approach, the AM is determined by the contribution of 
a single valence nucleon, proton or neutron. The only attempt to account for configuration 
mixing has been made in |7|, ||. However, as we shall show below, this is only the part of 
the many-body corrections, and is not the dominant one. The contribution arising from the 
induced PNC interaction in the nuclear core is considerably larger. It was not discussed yet 
at all. We present here the first treatment of these kind of effects. 

The major coherent effect induced by the residual interaction is the polarization of the 
nuclear core by a valence nucleon. This is the main effect causing the deviation of nuclear 
magnetic moments from the Schmidt values. The magnitude of polarization effects depends 
on the number of transitions from the core states over the Fermi surface. This number is 
determined by the selection rules, i.e. by the tensor rank of the operator. In the case of 
the AM, the number of transitions contributing to the core polarization is greater compared 
to the magnetic moment, and therefore, larger renormalization of the AM is expected. Just 
to illustrate the above statement, we refer to renormalization of the M3-octupole moment 
compared to Ml. In the case of an octupole moment, the number of transitions over the Fermi 
surface is much greater, and the core polarization reduces the valence nucleon contribution 
to M3 by a factor «4 | 

A convenient way to describe the core polarization is to use the effective renormalized 
operators, or effective fields, in the terminology of the theory of finite Fermi systems (1^ 



In the random phase approximation (RPA), the effective fields are the solutions of a system 
of integral equations describing particle-hole renormalization of the bare vertex. The weak 
nucleon-nucleon forces modify these equations. The modifications effectively produce an 
additional contribution to the AM compensating for the strong reduction of the single particle 
contribution. 

This paper is organized as follows. In the next two sections, we remind the set of opera- 
tors contributing to the AM, and discuss the accuracy of the leading approximation. Later 
on, we formulate the basic equations, and introduce the modifications of the equations by 
the weak nucleon-nucleon interaction in the leading approximation. Next, we calculate the 
renormalization of the single particle contribution, both analytically and numerically. Finally, 
we calculate all contributions to the AM both analytically and numerically, and discuss the 
stability of the results under variation of the constants of the strong residual nucleon-nucleon 
interaction. 
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2 The anapole moment operator 



The anapole moment operator is defined by |], ||, [TTJ 

a= -7r J d 3 r r 2 j(r), (1) 

where j(r) is the electromagnetic current density operator. 

The main contribution to the AM comes from the spin part of the current density. Nev- 
ertheless, the other contributions are noticeable and, apart from magnetization current, we 
shall discuss below the contributions from the convection, spin-orbit, and contact currents. 
Let us define the corresponding parts of the AM in the following way ||: 



a" = r x <j 

m 
lie 



m 



l ls 
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Here { , } is an anticommutator, p is the central nuclear density, f(r) = p(r)/p is the 
nuclear density profile, and Uf™ = U^f = 134 MeV ■ fm 5 is the proton-neutron constant of 
the effective spin-orbit residual interaction ||. The contact current contribution arises from 
velocity dependence of the effective nucleon-nucleon weak forces, taken in the form |5], [TT 



G 1 

Fw = "7h7" ({(9abcr a - g ba cr b ) ■ (p a - p 6 ), 5(r a - r b )} 

V^ 4m a ,b 

+ g' ab [* a x a b ] ■ V5(r a - r 6 )) = l - ^F w (ab). (3) 

ab 

Interaction (|3|) generates a mean field weak potential 

where # a = # ap ~ + ttn ~ . 

The effective interaction constants g ab , g ba , g' ab should be, strictly speaking, found from 
experiment. On the other hand, they can be estimated from the initial finite range PNC- 
interaction |T^] taking zero range limit with the account for short range particle particle 
repulsion [TB[. These "best values" estimates leads to g„ < 1, while the constant g p is 



approximately 4.5. The recent discussion of these constants [|14j give, however, different set 
for g p and g n , with g p ~ g n . Therefore, we shall keep below the constants explicitly as free 
parameters. 
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3 Single-particle contribution and 
leading approximation 

The leading approximation for the corrections to the single particle wave functions will be 
used in calculations of the core polarization effects. Therefore, it is worth to discuss the 
accuracy of this approximation. Neglecting the spin-orbit potential, and assuming constant 



nuclear density, we obtain for the correction to the single particle wave function [15] 

5ip (r) = -i£a(tr ■ r)^a(r), (5) 

where 

t - G 

Using (|^), we obtain for the spin part of the AM 

a ° = ^T— iRlrlR) WTTy (6) 

where K = (I — I) {21 + 1); R, I and I being the radial wave function, the orbital angular 
momentum of an outer nucleon and the nuclear spin. It is convenient to discuss AM in terms 
of a dimensionless constant k defined as (see Ref. 0) 

G Kl 

<ea>= 71^TT) K - (7) 

For the dimensionless constant k s we have 

k, = ^a^'\ (8) 

where we put {R\r 2 \R) = |r 2 ,4 2 / 3 . 

The naive expression (0) is in rather good agreement with exact numerical calculations 
of the spin part of AM [0. Therefore, the approximation is reasonable for averaging the 
volume type quantities like r 2 . For surface type quantities the situation is quite different. An 
instructive example is the convection current contribution to the AM. Numerical calculation 
gives for 209 Bi K conv = —0.019 while in the leading approximation (|5|) we obtain the value 
that differs by factor w 5. The explanation consists in the surface nature of the convection 
current contribution |J 

= -*g^{5R\ r 2 (-f + — ) \R)- (9) 

The integrand in the matrix element (^) for the outer nucleon with large angular momentum 
is peaked at the nuclear surface and the difference between exact 5R{r) and its approximate 
expression in leading approximation rR{r) provides considerable changes in the convection 
current contribution. 
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4 RPA renormalization of the AM 



The AM is a T-odd operator. Thus, the effective two particle interaction involved in AM 
renormalization must change sign under T-reversal of one of the two particles 

TaFiafyT- 1 = T b F{ab)T^ = -F(ab) 

The simplest interaction satisfying this condition is the same spin-spin interaction that 
changes nuclear magnetic moments: 

F s (ab) =C (g + g' T a ■ r b ) cr a ■ cr b 5(r a - r 6 ). (10) 

Here C is the normalization constant that we choose according to [|K| C = 300MeV ■ fm 3 
and the constants g' = 1.01 and g = 0.63. 

The effective interaction between the valence and core particles changes the interaction of 
the valence nucleon, with an external field producing additional core field. In the RPA this 
effect is accounted for by introducing a dressed effective vertex V satisfying the equation 

V = Vfcflai]) +FAV. (11) 

Here Vo([a^]) is one of the bare AM operators Eq.@; A is the static polarization loop of a 
particle-hole pair. 

/da 
— G UlI/2 {e)G u ' 2u/i {e), (12) 

G„ M (e) being a single particle nucleon propagator. In Eq.flTTD F is the sum of spin-spin 
interaction Eq.(|l0|) and the weak effective interaction Eq.(||). The propagator G VlV2 (e) should 
be calculated in the total mean field potential including the weak potential @ 

It is, however, more convenient to single out the weak interaction effects, treating them 
explicitly in first order perturbation theory. Let 5V be a correction to the vertex from the 
weak forces. For the unperturbed vertex V and the correction 5V we have the equations 

V = V Q ([a 4 ])+F 8 AV, (13) 

5V = F W AV + F s 5AV + F s A5V. (14) 
Here F w is the weak nucleon-nucleon interaction (|j). The AM value is given by 

a = (6il>\V\il>) + (il>\V\6il>) + (V>|<W)- (15) 

In the leading approximation, we have 

a = i£(iP\[*-r,V}\il;) + (iP\5V\iP). (16) 

The first term represents the single particle contribution renormalized by the spin-spin inter- 
action, while the second term is an additional contribution from the core particles. Note that 
the single-particle contribution is now the expectation value of a transformed renormalized 
anapole operator V. The AM is a T-odd operator of El type. The commutator of the AM 
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with i(j ■ r transforms it, as well as 5V, into a T-odd Ml type operator that evidently has 
nonzero expectation value in a state with spin /. The renormalization effects from the core 
polarization are different for these two types of operators. In the next two sections we shall 
discuss the renormalization of the electric type single particle AM operators and the magnetic 
type operators induced by PNC effects in the core. 

Note, that the contact term contribution produces the magnetic type operator from the 
very beginning. Therefore, its renormalization is similar to that of SV. 



5 Renormalization of the electric type single particle 
operator 

To solve Eq.fllBD, it is convenient to separate the angular dependence, introducing a set of 
tensor operators of rank J 

Tj M = {(T®Y Lm } JM . (17) 

The spin part (|2]) of the AM is proportional to T^ M . This is the only T-odd operator of 
the rank 1 with negative parity. Therefore, the dressed vertex V will have the same angular 
dependence as the bare vertex 

V s = v s {r)Tl M . 

The core polarization effects the radial dependence only, which for bare spin vertex is just 

vf(r) = N»r, 

where N" is 

Nf = 



' 87T 7ia[l a 



3 m 

The dressed vertexes v s (r) satisfy the following equations: 



u«(r) = vf{r) + V gf / r' 2 dr' A b {r, r')v b s (r'). (18) 

b=p,n J ° 

Here, the constants g v § and g^ 1 refer to the proton-proton and proton-neutron spin-spin 
interaction (|KD g^ = go + g' Q and g"^ = go — g' . The normalization interaction constant C 
is included in the radial polarization loop which in our case is 

A(r,r') = \c k 3ln \(jl\\Tl\\ 3 r)\ 2 R jln (r)R^ (19) 

jlnj'l' 

Here, kji n are the occupation numbers of filled nuclear states, {jl\\T±\\j'l') is the reduced 
matrix element of the tensor operator, and Gjqi{r, r'; tji n ) is the Green function of the radial 
Schrodinger equation with the angular momentum f, I' taken at the energy of the occupied 
level e jln . 
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Before going over to discussion of the numerical results, let us start from a simple model 
estimates of the core polarization in a harmonic oscillator potential without spin-orbit inter- 
action. In order to understand orders of magnitude, we shall calculate a polarization loop 
with the bare spin anapole vertex. Expanding symbolic Eq. flI3|) , we find for the first-order 
term 

a«(r) = CgJ^t(T)o'MT)*^('/\ — (* * (20) 

vv' u u 

Here, Vv(r) are the single-particle wave functions. In the absence of a spin-orbit potential 
one can sum over spin variables 

ai 1} W = -Cg^cr x £ ^ftM*) — ^(u'\r\u). 

vv u u 

In a harmonic oscillator we have a relation 

r=-— 2 [H,p}- (21) 



Using it, we obtain 



(i) iCg s vre/i 



a'/'lr) — — -tr x Y^ipl(r)i! u/ (r)(k u - k u/ )(u'\p\u) 

mar m , 

vv' 



muj z m muj A m 

Taking the expectation value of the correction in the state with total angular momentum I, 
we find for the ratio of this correction to the zero-order term 

< > _ < ^|(<T-r)aW|y; > 

< 5ip\a s \ip > < ip\(a- ■ r)sL s \ip > 

Cg sPo (R\rf'(r)\R) ^ 

muj 2 (R\r 2 \R) ~ ' 1 ' 

Since f'(r) is negative, the correction is negative and large. The sign of the correction is 
defined by the sign of the spin-spin interaction. The repulsive interaction decreases the single- 
particle contribution. For more realistic potentials, accounting for the spin-orbit potential, 
we can expect some changes in this ratio, since the correction is maximal on the nuclear 
surface, where the spin-orbit potential is important. Nevertheless, it remains large. 

The results of calculations of the renormalized single particle contribution are listed in 
Table I and Table II for proton and neutron levels. Note the reduction of the single particle 
contribution approximately by a factor of 2, in accordance with the above estimates. 

The spin-orbit and contact current contributions differ from the previous case only by the 
radial dependence of their bare vertex. Therefore, their renormalization can be done using 
the same Eq. (p~8|) . 



6 Renormalization of the magnetic type operators 



Let us now come back to Eq. (|14D describing the additional contribution to the AM coming 
from parity violation effects in the intermediate states of the core particles. Eq. ([TJ]) is of the 
same kind as (|13|) describing the renormalization of the valence nucleon contribution. The 
difference is in the driving force or the bare vertex, which is no longer connected to the bare 
AM operators (0). In calculation of the driving force, we shall use the leading approximation 

(!)■ 

Expanding the symbolic notation, we obtain 

F s (ab)5A b V b = dJ2 I d 3 r b ^t(r b )F s (ab)^(r b ) - - — (i/\V\u) 

up' J e »~ 6 »' 

= *6E / d\ b {^ u {v b )[ ( T b -v b ,F s {ab)}ij u ,{r b )(u'\V\u) 

uv' 

+ ^l{v b )F s {ab)^{nW\W ■ r, V\\u)\ ^— ^ 

> e„ — e u > 

= i£ b {[<T b ■ r b , F s (ab)} A b V b + F s {ab)A b [(T b ■ r b , V b ]} . (24) 

As we see, the driving force in Eq.(p4]) consists of two different parts. The first term can be 
combined with the weak interaction contribution F w (ab)A b V b , giving an effective interaction 
is the sum of the direct and inuced weak interactions 

F w (ab) = F w (ab) + z£ 6 [a b ■ r b , F s (ab)}. (25) 



The induced weak interaction was first introduced in fll6[ . It has the form 

F™ d {ab) = i[t a <T a ■ r a + i b a b ■ r b , F s (ab)), (26) 

and eventually appears in first order calculations of the residual iteraction. In our case, 
however, half of interaction (|26| ) enters equation fl24j) . Using the full interaction (p6|) will 
produce double counting, because the part related to the valence nucleon is already accounted 
for in the Sip correction to the valence nucleon wave function. The matrix elements of the 
induced weak interaction are proportional to the nuclear radius. Therefore, they are enhanced 
compared to the matrix elements of the direct weak interaction by the factor A 1 / 3 . For heavy 
nuclei, this is a considerable factor and, for this reason, we shall omit below the contribution 
of the direct term. 

The second term in ([24]) produces contributions of a different type to the AM. The com- 
mutator [er ■ r, V] produces the Ml type of operator, and we can expect its renormalization 
to be close to that of the magnetic moment. This very contribution has been previously 
discussed in J7|, ||. 

According to (p4|) , the correction SV can be presented as a sum of two terms satisfying 
but with different driving forces 

SV = 5V {1) + 5V {2 \ 
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5V^ = i£ b [<r b ■ r 6 , F s (ab)}A b V b + F s (ab)A b 5V b {1) 

5V^ = i£ b F s (ab)A b [(T b ■ r 6 , V b \ + F s (ab)A b 5V b {2) . (27) 
It is convenient to use, instead of 5V^ 2 \ another variable related to it via 

Xa = lU°a-Va,V a ]+5V^. (28) 

This variable satisfies the equation 

Xa = Oa • r a , V a ] + F s {ab)A bXb . (29) 

With this definition, we obtain from Eg. (|T6|) the following value for the anapole moment 

a= (il>\ X + 6VW\il>). (30) 

Thus, the contributions to the AM can be presented as the expectation value of the sum of 
two magnetic type operators induced by different driving forces. The first term represents 
the contribution of the renormalized magnetic operator obtained via Michel transformation 



15], while the second is the contribution of the induced weak interaction, similar to that of 



0. 

Let us now make an analytical estimate of this correction in the model used above for 
the estimates of the renormalization of single particle AM operator. We shall calculate the 
driving term in Eq.(^) for 8V^\ using instead of renormalized vertex V the bare spin vertex 
defined by Eq.(FJ). The correction to a proton contribution to the AM can be presented in 
the following form 

up _ up 

5a k ~ ^/vrJE^WK • r,o$M*W\(* * rr k )\u)^—^- 

vv' tv e v' 

+ ^f>n^E^(r)K • rX]W(r)<i/|(r x a^v)^— ^. (31) 



In Eq.(31), we have omitted factors common for protons and neutrons. Calculating the spin 
commutators we obtain 

5a k s ~ i P gfii p ai^{v) + ^f/^jf (r), (32) 

where 

A* k (r) = J>t(r)(r x <T)^ u ,(r)(v'\(r x V )*\ V )*LZ±L. (33) 

vv' e " €u ' 

Summing over the spin variables, we obtain for A lk (r) 

A ik (r) = 5 ik B ll (r)-B ik (r), 

where 

B ik (j) = ^l(T)r i MW\r k W)^^ L - (34) 

,,,,/ e v e v' 
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Using again relation ([H]), we obtain 



mar 



= — 2 Y,ti(?)[r\P k \M*)k v = l —5 ikP {v), (35) 

muj z „ muj z 

where p(r) is the proton or neutron density. Restoring the omitted factors, we have for the 
correction to the AM of a valence proton 

2neC 

5a s (r) = — {gf Hpipppir) + g p s n p n £ n p n (r)) cr (36) 

Its expectation value in a state with nuclear spin I is 

2neC K — - 

Sa s = — - (g PP Pp£pp P + gTPninpn) T(T ■ \X (37) 
m z uj z 1 [1 + 1J 



where K was defined in Eq.(||). For the "best values" of the weak interaction constants [12| 



in *C ip for heavy nuclei. Neglecting we obtain for the ratio of Sa s /a s for a valence proton 

6a, 2 Cg7Pa z fl _i^ 



a s muj z (R\r 2 \R) A V 2K 

where we use the total nuclear density po instead of the proton density. Using for u the 
standard value u = Al/A^-^MeV, we find 

5a s Z 

— «4.6— . 39 

a s A 

This result is quite instructive. The contribution to the AM from the induced weak interaction 
is greater than that coming from the single particle weak potential. However, one should keep 
in mind that the above calculation has been performed for the bare anapole vertex. 

7 Results 

The complete results of calculations are summarized in Table I and Table II. In the first 
column, we list the results of previous AM calculations in the independent particle model. 
In the second column, we list the contribution of the valence nucleon renormalized by the 
strong spin-spin interaction. On average, the core polarization reduces the single particle 
contribution by a factor w 2. In the third column, the AM induced by the core nucleons is 
written. It significantly restores the reduction of the single particle contribution, leaving the 
overall renormalization within pa 10%. 

The spin-spin residual interaction depends on the two constants go and g' corresponding 
to the interaction in isoscalar and isovector channels. The major part of the anapole moment 
is proportional to the nucleon magnetic moments; therefore, the isovector part of the AM 
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dominates. For this reason, we can expect small sensitivity of the AM to the isoscalar 
constant gr - F° r 205 T1 changing g in the interval 0.2 < g < 0.8 we get for the AM 0.313 < 
K to t < 0.380. The sensitivity to the isovector constant g' is larger, but still within reasonable 
limits. Changing g' in the interval 0.5 < g' Q < 1.5 we find for 205 T1 0.450 > n tot > 0.327. 
However, since the constant g' Q is fixed much better from the fit of magnetic moment, the 
above interval for K tot is in fact too large. 

To summarize, we have calculated the many body contributions to the nuclear anapole 
moment in the random phase approximation with effective nuclear forces. We found that 
the contribution to the AM from the valence nucleon is reduced by a factor ps 2 from the 
core polarization effects. However, PNC effects in the core states produce an additional 
contribution to the AM that partially compensates the reduction of the single particle AM. 
The resulting value of the AM appears to be close to its initial single-particle value calculated 
in the independent particle model. 
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Table 1: Anapole moment contributions for proton levels 
s.p. (5i/j\V\i/j) WW) 



Total 



K s 

Kls 

K-conv 
K c 

Ktot 



0.070< 
-0.017<£ 
-0.0041^ 
0.0066^; 



pit 

w 



0.050< 
-0.0066<™ 



0.039g£ 
-0.0094< 
-0.0062< 



0.024< 



0.018^+0.009^ 
-0.0035^-0.0004^ 
0.0037^-0.0001^ 
0.0052<™-0.0006# 



np 



0.018^+0.008^ 
+0.0052^-0.0006^ 



w 
np 



0.057 <+0.009< 

- 0.013^-0.0004^ 
- 0.0025^-0.0001^ 
0.0052 <™-0.0006#™ p 

0.041 <+0.008< 
+0.0052<™-0.0006c/™ p 



Kls 
K r 



K tot 



0.108^ 
-0.022<& 

-o.on< 

0.0085<,™ 

0.075< 
+0.0085^" 



0.052<£ 
-0.01L& 
-0.009< 



0.032^ 



0.048^+0.001^ 

-0.009^-0.0001^ 
0.0009^+0.00002^ 
0.0064<,™-0.0006#™ p 



0.039<- 
-0.0064<,™- 



-0.0003^ 
0.0006#™ p 



0.100 ^+0.001^ 
- 0.020^-0.0001^ 
-0.0077^+0.00002^ 
0.0064 #£, n -0.0006#™ p 

-0.018 <-0.0003#™ 
+0.0064<,™-0.0006c/™ p 



^•conv 
K c 

Ktot 



0.083< 
-0.024^ 

-0.005< 

o.oio<c 



0.039< 
-0.011 ff £ 
0.0009< 



0.032<£. 
-0.007<& 
0.003^ 
0.006<,™- 



f0.0045#™ 
-0.0003<; 
0.00004#™ 
-0.0004^ 



0.055< 
+0.010<" 



0.029# 



0.028^+0.004^ 
+0.006<™-0.0004< p 



0.071^+0.0045^ 
-0.020^-0.0003^ 
0.004^+0.00004^ 
0.006^ n -0.0004^ 

0.057 <+0.004< 
+0.006<™-0.0004< p 



Table 2: Anapole moment contributions for neutron levels 







s.p. 


(5i>\V\i>) 




(i>\6V\i>) 


Total 




Kg 


-0.099< 


-0.052,?™ 


-0.0035< 


-0.028(7™ 


-0.0035^-0.080^™ 




Kls 


0.006^ 


0.003^ 


0.0008< 


+0.002(7™ 


0.0008 <+0.005<, 


207p b 


^conv 




-0.0002^ 


0.0001<&- 


■0.0002(7™ 


0.0001^-0.0004^™ 




K c 


-0.002^ 




0.0001^™- 


-0.001<f 


0.0001 <™-0.001#™ p 




Ktot 


-0.093(7™ 


0.050(?™ 


0.003<?£ 


-0.027(7™ 


-0.003 <-0.076(7™ 






-0.002<f 




+0.0001^™- 


-0.001^ 


+0.0001<,™-0.001<f 




Kg 


-0.062^ 


-0.033,?™ 


-0.004< 


-0.033(7™ 


-0.004 <-0.066(7™ 




Kls 


0.01< 


0.004,?™ 


0.001<+0.002< 


0.001 <+0.006(?™ 


209p b 


^•conv 




-0.0003(7™ 


-0.0005<,- 


-0.0002(7™ 


-0.0005^-0.0006^™ 




K c 


0.004^ 




0.0003^™- 


-0.003^ 


0.0003 ^™-0.003c/™ p 




Ktot 


-0.052,?™ 


-0.029(7™ 


-0.004< 


-0.031(7™ 


-0.004 <-0.060< 






0.004^ 




+0.0003<™- 


-0.003<? 


+0.0003^-0.003 gl p 
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